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Abstract. We provide probabilistic interpretation of resonant states. This we do by showing that 
the integral of the modulus square of resonance wave functions (i.e., the conventional norm) over a 
properly expanding spatial domain is independent of time, and therefore leads to probability conser- 
vation. This is in contrast with the conventional employment of a bi-orthogonal basis that precludes 
probabilistic interpretation, since wave functions of resonant states diverge exponentially in space. 
On the other hand, resonant states decay exponentially in time, because momentum leaks out of the 
central scattering area. This momentum leakage is also the reason for the spatial exponential diver- 
gence of resonant state. It is by combining the opposite temporal and spatial behaviors of resonant 
states that we arrive at our probabilistic interpretation of these states. The physical need to normalize 
resonant wave functions over an expanding spatial domain arises because particles leak out of the 
region which contains the potential range and escape to infinity, and one has to include them in the 
total count of particle number. 
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1. Introduction 

The present paper reviews the basic definition of the resonant state in quantum mechan- 
ics [1] as well as a recent progress on physical interpretation of the resonant state [2]. 
Presenting partly new calculations, we stress that the particle number, after its proper def- 
inition, is conserved for the resonant state. Hence, the probabilistic interpretation of the 
square modulus of the wave function remains intact even for the resonant state. 

The phenomenon of resonance has been of great importance in elementary-particle 
physics and nuclear physics [3-28]. In recent years, unstable nuclei, which are nothing 
but resonant states, are experimentally sought after; see e.g., Refs. [29,30]. It has been also 
realized in condensed-matter physics that the resonance reveals itself in electronic conduc- 
tion in mesoscopic systems [31-36]. These facts, among others, have revived theoretical 
interest in the phenomenon of resonance in quantum mechanics [37]. 

The wave function of a resonant state has been discussed greatly since the early stage 
of quantum mechanics. As will be shown at the end of Sec. 3 below, it exponentially di- 
verges in the distance from the scattering center. Although this is an inevitable conclusion, 
the divergence has hindered the probabilistic interpretation of the square modulus of the 
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wave function. A bi-orthonormal space of the left- and right-eigenfunction was thereby 
introduced to define a convergent norm [7,38^-5]. 

The main purpose of the present paper is to stress that the probabilistic interpretation of 
the square modulus of the wave function is nonetheless possible [2]. We show that particles 
leak from the central region in the resonant state. The exponential divergence in space of 
the wave function is a direct consequence of the leaking particles. The particle number (in 
the form of the square modulus of the wave function) is conserved when we take the leak 
into account. More specifically, we expand the region where we count the particle number, 
with the speed of the leaking particles. For arbitrary superpositions of states, we introduce 
a definition of the speed of the leaking particles. With this proper definition of the speed, 
we show that the integral of the square modulus of the wave function over the expanding 
region is independent of time. 

The paper is organized as follows. In Sec. 2, we review the definition of the resonant 
state as an eigenstate of the Schrodinger equation. We here introduce the Siegert condition 
that there exist only out-going waves away from the scattering region. We then present 
in Sec. 3 a physical view of the resonant state with the Siegert condition. The resonant 
state will be described as a state that decays in time exponentially with momentum leaks 
from the scattering region. On the basis of this view, we show in Sec. 4 that the particle 
number defined as the square modulus of the wave function is indeed conserved. The 
probabilistic interpretation of the square modulus of the wave function is thus intact even 
for the resonant state. The last section is devoted to a summary. 

2. Direct definition of quantum resonance 

In many textbooks, resonance is defined as a pole of the S matrix. It is, however, more 
convenient to define a resonant state as an eigenstate of the Schrodinger equation with the 
Siegert condition given below. Let us briefly review the connection of the two definitions 
in the present section. 

Consider, for simplicity, a one-dimensional system with a scattering potential localized 
around the origin. An incident wave Ae lkx from the left (Re k > 0) will be scattered to 
result in a reflection wave Be~ lkx and a transmission wave Ce lkx . When we properly 
connect these three waves at the scattering potential, the coefficients A, B and C become 
dependent on the wave number k, or on the energy E through the dispersion relation E(k) 
of a free space (specifically E{k) = h 2 k 2 / (2m) in the following example). 

Elements of the S matrix are then given by the ratios B/A and C/A. The resonance 
is defined as a pole of the S matrix in the complex wave number plane or in the complex 
energy plane. In fact, zeros of A give poles of the S matrix. 

The above two-step procedure of (i) solving a scattering problem for real k and E and 
(ii) finding zeros of A in the complex k and E plane, can be reduced to just one step of 
solving a scattering problem for complex k and E with the condition A = required from 
the very beginning. The condition A = means that we have out-going waves only. This 
boundary condition is referred to as the Siegert condition. The above argument shows that 
a pole of the S matrix is indeed given by a state with the Siegert boundary condition of 
out-going waves only [4,5]. 

For some readers who might wonder if a solution with the Siegert condition is possible, 
let us solve a simple problem shown in Fig. 1. The Hamiltonian is given by 
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Figure 1. A square-well potential with the Siegert condition of out-going waves only. 



with the square-well potential 

^) = {„-"» <0 

We assume the form of the solution as 

Be~ %Kx forx<-2, 
$(x) = { Fe~ lK ' x + Ge lK ' x for \x\ < I, (3) 
Ce lKx for x > I, 

where K and K' are wave numbers in the respective region, related to each other as 

- 2 ™V = K\ (4) 
n 

We here used uppercase letters for the wave numbers in order to emphasize that they are 
complex numbers. 

Since the potential is an even function, we can classify resonant solutions into odd and 
even ones. We set B = —C and F = — G for odd solutions and B = C and F = G for 
even solutions. Upon solving the Schrodinger equation, we require the usual conditions 
that $>(x) and $>'(x) are continuous at x = I. First, the continuity of at x = I gives 
odd solutions of the form 

® [x) - 2tGX {sgnx S m(K>l)e iK ^- 1 ) for \x\ > I, (5) 
and even solutions of the form 

*{x) - 2G x | cos{K , iy K(\^i) for ^1 > i (6) 

Then, the continuity of ^'(x) at |x| = / gives the equations 
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Figure 2. Solutions K n and E n for the square-well potential with the Siegert condi- 
tion. We set Vb = l{h 2 /2ml 2 ]. The dot on the positive imaginary K axis and on the 
negative real E axis is a bound state. (It has a small real negative value of the energy.) 
The dots on the fourth quadrant of the K plane and on the lower half of (the second 
Riemann sheet of) the E plane are resonant states, whereas the dots on the third quad- 
rant of the K plane and on the upper half of (the second Riemann sheet of) the E plane 
are anti-resonant states. 



K = -iK' cot(K'l) for odd solutions, (7) 
K = iK' tan(K'l) for even solutions, (8) 

which should be solved together with Eq. (4). We solved them numerically to obtain K n 
and E n plotted in Fig. 2. First of all, the wave function of the form (3) represent a bound 
state when K is a positive pure imaginary. We then have pairs of solutions on the lower 
half of the K plane. One of each pair on the fourth quadrant is called a resonant state, 
whereas the other one on the third quadrant is called an anti-resonant state. (In the context 
of electronic conduction of mesoscopic systems, some call a dip in the conductance an anti- 
resonance. This is completely unrelated to the anti-resonant state in the present paper.) 

Some readers might wonder why the present potential without any potential barriers 
can support resonances. Indeed, a conventional view of quantum-mechanical resonance 
might be as follows for a potential with some barriers such as exemplified in Fig. 3(a); a 
particle that would be completely trapped classical-mechanically in the range \x\ < I and 
hence would be a bound state quantum-mechanically, can in fact escape the trap because of 
tunneling effect, and this coupling to the continuum outside the barriers causes a resonance 
width to the would-be bound state; see Ref. [2] for such an example. Such a resonance 
shows up in the Hermitian problem as narrow transmission resonances at energies below 
the barrier, where the transmission coefficient shoots up to unity. 

Nevertheless, the potential well in Fig. 3(b) can support resonances that we discuss 
here and hereafter. The square-well potential can partially trap a quantum-mechanical 
particle because the wave of a particle that departs the origin x = towards outside is 
partially reflected at the potential edges x — ±1 even when its kinetic energy is posi- 
tive. This reflection at the potential edges plays a role of an effective potential barrier. 
Indeed, this effect would cause the Anderson localization if we had infinitely many po- 
tential wells [46]. We will argue in the next section that the complex eigen-wave-number 
and eigenenergy of a resonant state are directly linked to out-going momentum fluxes from 
the central range of the scattering potential. This view holds for potential of both types 
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Figure 3. (a) A particle trapped in the shaded area can escape the trap because of tun- 
neling effect. The coupling between the bound state inside the trap and the continuum 
outside makes the bound state unstable and results in a resonance width, (b) A particle 
with a positive kinetic energy is partially reflected at the potential edges x = ±Z and 
hence can be "trapped" in the potential well for a while. 



in Figs. 3(a) and (b). The transmission coefficient associated with the square-well poten- 
tial in Fig. 3(b) is very close to unity for any positive energies with broad local maxima 
around where it reaches unity. These are not well-defined transmission resonances in the 
conventional sense. Nonetheless, they are in fact the Hermitian shadows of the resonances 
that we discuss in the present paper. Interested readers might want to follow the discussion 
hereafter with their own example potential of the type in Fig. 3(a); it will not cause much 
difference. 

At the end of the present section, we stress that the above definition on the basis of 
the Siegert condition gives the resonant (and anti-resonant) states as eigensolutions of the 
Schrodinger equation. The wave function of each resonant or anti-resonant state is explic- 
itly given by substituting the corresponding solution in the K plane in Eq. (5) or (6). 



3. Physical view of quantum resonance 

In the present section, we give a physical view of the resonant and anti-resonant states 
as defined in the previous section. We begin with discussing the non-Hermiticity of 
the Hamiltonian of open quantum systems. For simplicity again, we consider a one- 
dimensional system of the form (1), where V{x) is a real function with a compact support 
around the origin. 

We first define the expectation value of the Hamiltonian with respect to an arbitrary 
function ^>(x) in the form 

(*|W|*)n = J V{x)*m>{x)dx, (9) 

where £1 = [-L, L] is a large region whose edges x = ±L are far away from the support 
of the scattering potential V(x); that is, V(L) = V(—L) = 0. We may take the limit 
L — ► oo in the end if we can, but for the moment we keep L to be finite. (Note here and 
hereafter that the bra vector is given by the standard definition of being complex conjugate, 
not given in the dual space spanned by the left eigenvectors [7,38-45].) 
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We are now interested in the non-Hermiticity of the Hamiltonian; that is, we would like 
to know whether (^E'|H|5')o is real or not. Clearly, (\D , |V(x)|\f r )£2 is real. The remaining 
problem is the reality of the quantity (^\(d 2 /dx 2 )\^)n. For this purpose, we carry out a 
partial integration: 



dx 2 



^ = J *(x)**"(x)dx 



= [f(x)*M>'(x)]^ = _ L - f \V(x)\ 2 dx. (10) 

The second term on the right-hand side is obviously real. On the other hand, we rewrite 
the first term as follows; 

where p = (H/i)d/dx is the momentum operator. Therefore, we have the identity 

Im<*|W|*>n = -^Re{^\p n \^} dn , (12) 

where we used the bra-ket notation, which holds for any space dimension d. Here, p n is 
the normal component of the momentum on the boundary of 0; in the present case, p n = p 
at x = L and p n = — p at x = -L. We also let (*| • \^)dn denote the expectation value on 
the boundary of 0. We derived the identity (12) in one dimension but it holds in arbitrary 
dimensions, in fact. (When we put $ to each resonant state $„, the identity (12) reduces 
to the one derived previously [47-49].) 

The equality (12) shows that the non-Hermiticity of the Hamiltonian comes from the 
momentum leak at the boundary of the region. Its meaning becomes clearer when we ex- 
tend the above argument to solutions of the time-dependent Schrodinger equation, *&(x, t). 
For a solution of the equation 

ih^*(x,t) = H*(x,t), (13) 
we can easily prove the identity [50] 

^<*|*> n = |lm(*|W|*)n. (14) 
and thus we have 

dt m 

This has a very plausible interpretation; the decreasing rate (for Re('J'|p n |^')an > 0) of 
the particle numbers in the region fl is proportional to the momentum leak at the boundary 
of the region, d£l. Equations (12) and (14) show that the two quantities are indeed linked 
to the non-Hermiticity of the system. This is how an open quantum system becomes non- 
Hermitian. 
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Figure 4. A schematic view of the particle number conservation of a resonant state. 
The particle number decays in time as exp(— Im E n t/Ti) over the entire space, but it 
is conserved when we expand the integration region with the speed v n — Re K n h/m. 



The above interpretation leads to the following physical view of the resonant and anti- 
resonant states. In a resonant state, there is a leak on the boundary of the region ft 
(ReK n > 0) and hence the particle number in il decays in time (ImE n < 0). In an 
anti-resonant state, on the other hand, particles come into the region £1 (Re K n < 0) and 
hence the particle number in il grows in time (Im E n > 0). This view of the resonant and 
anti-resonant states applies both to a trapping potential in Fig. 3(a) and to a potential in 
Fig. 3(b). 

Incidentally, the dispersion relation E n = (HK n ) 2 /2m is followed by 

h 2 

ImE n = — ReK n 1mK n . (16) 
m 

The argument in the previous paragraph then gives Im K n < for all resonant and anti- 
resonant states. This, however, means that the wave function (3) diverges in the form 

\$ n (x)\ 2 „ e 2|iw*»IM as|a;|-»oo. (17) 

This divergence indeed has hindered the probabilistic interpretation of the wave function 
of the resonant state. Nevertheless, we show in the next section that the probabilistic 
interpretation is actually possible. 



4. Particle number conservation 



In the present section, we show that the particle number is conserved for the resonant and 
anti-resonant states. For simplicity again, we proceed in one dimension. The basic idea is 
as follows (Fig. 4). Consider a resonant eigenstate 

$ n (x,t) =$ n {x)e- iE ^ h (18) 

with the eigen-wave-number K n and the eigenenergy E n . At one point of time, which we 
refer to as t = 0, we count the particle number of the state over a finite region of space, 
fl = [— L(0), L(0)]. In order to compare the particle number at a later time, we have to 



Pramana - J. Phys., Vol. xx, No. x, xxx xxxx 



7 



N Hatano, T Kawamoto and J Feinberg 



take into account the fact that the particles are fleeing from the region in a resonant state, 
as we pointed out at the end of the previous section. The fleeing speed outside the potential 
range may be considered as 

v n = -ReK n . (19) 

m 

Therefore, we should expand the region Cl(t) = [— L(t), L(t)] for counting the particle 
number. The speed of the expansion should be equal to v n ; that is, 

L(t)=v n t + L{0). (20) 

We show below that the particle number in the expanding region is indeed constant [2]. 
Let us calculate the time derivative of the particle number in the expanding region: 




= -^Re<*|A,|*> 8 n(t) +2L(t) \^(L(t),t)\ 2 

= ~Re{*\p\*) x=L(t) + 2L(t)<*|*) x=i(t) . (21) 

(Note that we do not take space integration in both terms of the last line. In higher di- 
mensions, we would need surface integration over the boundary of the region Q.) The first 
term on the right-hand side of Eq. (21) is given by Eq. (15) and is equal to the decay rate 
of the particle number in the region il, while the second term comes from the derivative of 
the integration bounds ±L(£) and is equal to the increase of the particle number due to the 
expanding integration region. 

Using the explicit form of the wave function (3), or 

! e i(K n x-E n t/h) 
on the right of the potential range, 
±e i { -K nX -E n t/h) ( 22 ) 
on the left of the potential range, 

where the sign of the second line on the right-hand side depends on whether the solution is 
odd or even, we obtain 

j t (*n|*n> n(t) = 2 (-^RcX„ + L{t)j e M*»Kt)-*„t/ft) = o. (23) 

Not only the coefficient but also the time dependent part of the exponent on the right-hand 
side vanishes because of the dispersion relation E n — (hK n ) 2 /(2m). Thus we proved 
that each of the resonant and anti-resonant states conserves the particle number. In other 
words, the probabilistic interpretation of the square modulus of the wave function is intact 
for the resonant and anti-resonant states. 

The probabilistic interpretation of pure resonant states should also hold for superposi- 
tions of resonant and anti-resonant states. Thus, let ^(x, t) be a generic solution of the 
time-dependent Schrodinger equation (13). Since the identity (21) holds for arbitrary such 
functions <£, we can make (d/dt) (^ , |^ , )o( t ) vanish by requiring L(t) to satisfy 
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L(t) = v(L(t)) (24) 
with the speed of the fleeing particles defined by 

_ i (25) 

which generalizes Eq. (19). Equation (24) constitutes a differential equation for L(t). 
Thus, given t), we can determine L(i) by solving Eq. (24) instead of using Eq. (20). 

As a simple application of Eq. (25), beyond the pure resonant case, consider a superpo- 
sition of two resonant states 

^{x,t) = oi$i(o;,t) + a 2 $ 2 (x,t) (26) 

with time-independent coefficients a\ and a 2 . Outside the potential range, we have the 
form (22) for each of $1 and 3> 2 with ReK n > 0, lmK n < and lmE n < 0. Substi- 
tuting of Eq. (26) in Eq. (25) then yields, after a straightforward (but somewhat tedious) 
computation, to the differential equation 

• h/m 



| ai |2 e ImA(t) + | a2 |2 e -ImA(t) + CQsQ ^ 

x {|ai| 2 e ImA « Retfi + \a 2 \ 2 e- ImA ^R C K 2 
+ \ ai a 2 \ [ReiKx +K 2 )cose(t) +lm(K 1 - if 2 ) sin6(t)]}, (27) 



where 



A(t) = El - E2 t - (Kj. - K 2 )L{t), (28) 
Q(t) = ReA(t) - (argai - arga 2 ). (29) 

Let us solve the differential equation (27) for the first and second resonant states of the 
symmetric square well potential (2). The states are depicted as the first and second dots in 
the fourth quadrant of Fig. 2(a). Specifically, we have 

K x ~ 2.356987 - il.909078[l/Z] 

with E 1 ~ 1.910812 -i8.999349[?j 2 / 2 wZ 2 ], (30) 
K 2 ~ 4.119962 - i2.301222[l/Z] 

with E 2 ~ 11.678469 - il8.961903[?i 2 /2m/ 2 ] (31) 

for Vq = l[K /2ml 2 ]. We also use ai = a 2 = 1 for simplicity. The fleeing speed (25) 
depends on x as shown in Fig. 5(a). The speed is negative near the origin, but is posi- 
tive all through the region \x\ > I, converging to Re K 2 h/m fairly quickly. We show in 
Fig. 5(b) the numerical solution L(t) and L(t) of the differential equation (27) with the 
initial condition L(0) = I. The solution L{t) behaves almost linearly after some time with 
L(t) converging to Re K 2 K/m. This demonstrates that the resonance with greater | Im K n 
quickly dominates the behavior. 
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Figure 5. (a) Enumeration of the fleeing speed (25). The shaded area indicates the 
potential well, (b) The numerical solution of the differential equation (27) with the 
initial condition L(0) = I. 



5. Summary 

In the present paper, we first reviewed the definition of the resonant and anti -resonant states 
as eigenstates of the time-independent Schrodinger equation. These states are eigenstates 
with the Siegert condition, that is, the boundary condition that there are out-going waves 
only. We demonstrated for a square-well potential that the resonant and anti-resonant states 
are indeed obtained under the Siegert condition. 

We then presented a physical view of the resonant and anti-resonant states. In the res- 
onant states, the particle number around the central region exponentially decays in time 
because of momentum leaks from the region. The anti-resonant states are the time re- 
versal of the resonant states. In this sense, the pair of a resonant and anti-resonant states 
spontaneously breaks the time-reversal symmetry [51-53]. 

The above physical view led us to the argument of the particle-number conservation 
for the resonant and anti-resonant states. In either case, one must consider both spatial 
and temporal behavior of the resonance wave function in order to have a probabilistic 
interpretation thereof. When we expand the integration region following the momentum 
leaks, the particle number in the expanding region is conserved. We then extended the 
probabilistic interpretation to a superposition of resonances. The latter extension, as well 
as details of the calculation in Sec. 4 appear here for the first time. Finally, we comment that 
the next obvious step would be have a clear probabilistic interpretation of superpositions 
of bound states and resonant states. 

The work of NH is supported by Grant-in-Aid for Scientific Research No. 17340115 
from the Ministry of Education, Culture, Sports, Science and Technology as well as by 
Core Research for Evolutional Science and Technology (CREST) of Japan Science and 
Technology Agency. The work of JF was supported in part by the Israel Science Founda- 
tion (ISF). 
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